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ABSTRACT 

V921 Scorpii is a close binary system (separation 0.025") showing the B[e] -phenomenon. The system 
is surrounded by an enigmatic bipolar nebula, which might have been shaped by episodic mass-loss 
event s, possibly triggere d by dynamical interactions between the companion and the circumprimary 
disk JK raus et aT1l2012a[) . In this paper, we investigate the spatial structure and kinematics of the cir- 
cumprimary disk, with the aim to obtain new insights into the still strongly debated evolutionary stage. 
For this purpose, we combine, for the first time, infrared spectro-interfcromctry (VLTI/ AMBER, 
A/AA = 12,000) and spectro-astrometry (VLT/CRIRES, A/AA = 100,000), which allows us to study 
the AU-scale distribution of circumstellar gas and dust with an unprecedented velocity resolution 
of 3 kms . Using a model-independent photocenter analysis technique, we find that the Br7-line 
emitting gas rotates in the same plane as the dust disk. We can reproduce the wavelength-differential 
visibilities and phases and the double-peaked line profile using a Keplerian-rotating disk model. The 
derived mass of the central star is 5.4±0.4 M© • (d/1150pc), which is considerably lower than expected 
from the spectral classification, suggesting that V921 Sco might be more distant (d ~ 2 kpc) than com- 
monly assumed. Using the geometric information provided by our Br7 spcctro-intcrferometric data 
and Paschcn, Brackett, and Pfund line decrement measurements in 61 hydrogen recombination line 
transitions, we derive the density of the line-emitting gas (N e = 2. ..6 x 10 19 m -3 ). Given that our mea- 
surements can be reproduced with a Kcplcrian velocity field without outflowing velocity component 
and the non-detection of age-indicating spectroscopic diagnostics, our study provides new evidence 
for the pre-main-sequence nature of V921 Sco. 

Subject headings: stars: pre-main sequence — stars: individual (V921 Scorpii) — stars: emission- 
line, Be — binaries: close — protoplanetary disks — accretion, accretion disks — 
techniques: interferometric 



1. INTRODUCTION 

The class of the B[c] stars comprises some of the 
most peculiar and mysterious objects of stellar astro- 
physics. In particular, it has been found that the class- 
defining characteristics, which include infrared excess 
and strong permitted and forbidden line emission, are 
observed in a wide range of evolutionary stages, includ- 
ing stars of pre-main-sequence, post-main- sequence, and 
unkn own nature (unclassified B[e] stars, ILamers et al.1 
|1998[) . An enigmatic member of the latter category is 
V921 Scorpii (=Hen 3-1300, MWC865, CD-42°11721). 
Despite more than one hundred studies that have been 
conducted over the last four decades, the stellar pa- 
rameters and evolutionary stage of V921Sco remain 
strongly debated. These uncertainties arise primar- 
ily from problem s in distance estimation , where values 
between 160 pc (Hillenbrand et all I1992D and 2.6 kpc 
(|Shore et al.l I1990D have been proposed. Furthermore, 
due to the lack of photospheric lines in the opti- 
cal spectrum, the existing spectral classifications are 
based on rather indirect methods, such as photomet- 
ric modeling or the ionization analysis of circumstel- 
lar emission lines. Accordingly, the derived effective 

1 Based on observations made with ESO telescopes at the 
Paranal Observatory under the open-time program ID 084. C- 
0668(A, B) and with the Magellan Clay telescope. 



temp e ratures (12,300 K to 31,600 K; iHillenbrand et al.1 
Il992t iCidale et all 1200 ID and s pectral types (Ae p to 
B0[e]p: Ide Winter fe Thel [19901 iNatta et al.1 Il99l ex- 
hibit a considerable spread, which prevents a clear 
Hertzsprung-Russell diagram classification and has lead 
authors to argue both for an evolved (supergiant; 
Hutsemekers fc van Dromlll990l : Borges Fernandes et al.1 



20071) and youn g fHerbig B[el: fdTWinter fe Thd 119901: 
The et al.l 119941: iBenedettini et al.l 119981: lHabart et al.1 



2003tlAcke eTaLll2005l : IAcke fc van den Anckerll2006[ ) ob- 
ject nature. 

The aim of this paper is to gain new insights on the 
astrophysical nature of V921Sco, and B[e] stars in gen- 
eral, by studying the spatial distribution and kinemat- 
ics of the circumstellar gas and dust. These constraints 
will inform us about the disk-formation mechanism (i.e. 
accretion versus excretion), and, thus, the evolution- 
ary state of V921 Sco. For instance, accretion disks 
around young stars are believed to exhibit a Keplerian 
rotation profile, while the radiation-driven winds from 
evolved stars s hould exhibit a strong outf lowing veloc- 
ity component (jLamers fc PauldrachUL99l . Given the 
possible kiloparsec distance of V921 Sco, this task re- 
quires both high angular and high spectral resolution, 
which we achieve in our study by combining, for the 
first time, two highly complementary techniques, namely 
spectro-interferometry (VLTI/ AMBER, providing mil- 
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Fig. 1. — tre-coverage of our VLTI/ AMBER observations us- 
ing the LR-HK mode (black dots, Paper I) and HR-K mode (red 
squared, this paper), as well as the slit orientation of the spectro- 
astromctric observations (dashed grey lines). 

liarcsecond angular resolution and spectral resolution 
R = 12,000) and spectro-astrometry (VLT/CRIRES, 
R = 100,000). From these spatially and spectrally re- 
solved constraints, we derive the gas velocity field on 
scales down to a few stellar radii . 

In a recent study (jKraus et alll2012al Paper I), we al- 
ready obtained milliarcsecond-resolution continuum in- 
terfcrometric images of V921 Sco in three wavelength 
bands (1.65, 2.0, and 2.3 /jm) and discovered a close (sep- 
aration p ~ 25 milliarcseconds=mas) companion, as well 
as a continuum-emitting disk with an apparent size of 
7.5±0.2 mas (Gaussian FWHM along major axis). In ad- 
dition, we obtained images of the surrounding large-scale 
bipolar nebula and detected multi-layered, shell-like sub- 
structures that might have been shaped by episodic mass- 
loss events. Based on roughly matching timescales be- 
tween the estimated orbital period and the mass-ejection 
period, we suggested that the mass-ejection events might 
be triggered by the newly discovered companion. 

For the stellar parameters of V921 Sco, we adopt in 
this study a luminosity = (10 ± 3) x 10 3 L Q , effective 
temperature T c ft = (14 ± 1) x 10 3 K, a distance of d = 
1150 ± 150 pc, Ay = 4.8 ± 0.2, and a stellar rad ius of 
= 17.3 ±0.6 R Q ijBorges Fernandes et al.|[2007h . The 
implications of a modified set of fundamental parameters 
will be discussed in Sect. 15.21 

In the following, we present our spectro-intcrfcrometric 
(Sect. I2.1j) . spe ctro -astrometric (Sect. I2.2j) . and spec- 
troscopic (Sect. 12. 3p observations, followed by a discus- 
sion and quantitative modeling, both in the continuum 
(Sect. [3]) and the Br7-line (Sect. [4]). Finally, we will in- 
terpret our modeling results (Sect. [5]) and conclude with 
a brief summary of our findings (Sect. |6]). 



2.1. 



2. OBSERVATIONS 
VLTI/ AMBER spectro-interferometry 



Infrared interferometers, such as ESO's Very Large 
Telescope Interferometer (VLTI), coherently combine the 
light from separate apertures in order to achieve an un- 
precedented angular resolution of a few milliarcseconds 
at infrared wavelengths. With a spectral resolving power 
of up to R — A/AA = 12,000, these instruments now 
enable investigations both in the continuum emission 
and in spatially and spectrally resolved gas-tracing lines. 
The measured visibility amplitudes, closure phases, and 
wavelength-differential phases (DPs) provide powerful 
constraints for model-fitting or can be used for the re- 
construction of interferometric images. 

First spectro-interferometric observations on V921 Sco 
using the VLTI near-infrared bea m combiner instru- 
ment AMBER (iPetrqv et all 120071) were presented by 
iKraus et all (|2008bf ). iKreplin et all (|2012f ). and in Pa- 
per I and provided low (LR-HK, R = 35) and medium 
(MR-K, R = 1500) spectral resolution. Here, we present 
new observations using AMBER's unique high spectral 
resolution mode (HR-K, R = 12,000). The new data 
sets were obtained on 2010-03-03 and 2010-08-26 using 
the 8.2 m unit telescopes UT2, UT3, and UT4 and cover 
a wavelength window around the Br7 2.16602 /im line. 

In order to yield sufficient signal-to-noise ratio (SNR) 
for our high spectral resolution observation s, we em- 
ployed the FINITO fringe trac ker instrument (|Gai et all 
12004: iLe Bouauin et all [20081) . This instrument mea- 
sures and corrects the atmosphere-induced phase pertur- 
bations and allowed us to record the optical path delay 
(OPD)-stabilizcd interferograms with a long detector in- 
tegration time of 3 s. Unfortunately, the tracking per- 
formance of FINITO was degraded for the longest base- 
lines due to the low visibility contrast in this baseline 
regime. Based on the poor SNR, we rejected the vis- 
ibility measurements for one baseline from 2010-03-03 
(UT2-UT4) and two baselines from 2010-08-26 (UT2- 
UT4, UT3-UT4). The remaining baselines exhibit fringe 
SNRs up to 3.5, 4.5, and 14, and cover baseline lengths 
between 36.7 and 84.3 m and PAs between 22 and 292° 
(Fig. [IJ. The data shows a double-peaked, rising vis- 
ibility profile in the Br7-line as well as non-zero DPs 
(Fig. EJ . 

The wavelength calibration was done using atmo- 
spheric telluric features close to the B^-line and by ap- 
plying a heliocentric-barycentric system correction with 
heliocentric velocities of +28.33 kms" 1 (2010-03-03) and 
-27.09 km" 1 (2010-08-26), respectively. 

Each observation on V921 Sco was accompanied by ob- 
servations on the interferometric calibrator HD 161068, 
for which we assume a uniform -disk diameter of 1.448 ± 
0.019 mas (Merand et al 2006). Both for the science and 
calibrator star observations, we extract raw visibilities, 
closure phases, a nd DPs using the amdlib (V3.0) data re - 
duction software (jTatuMi et alll2007HlChelli et al.ll2009ft . 

Due to the influence of residual vibration, it is known 
that AMBER+FINITO observations with long DIT, 
such as used for our HR-K mode observations, might 
result in a poor absolute visiblity calibration, while the 
wavelength-di fferential observable s are only marginally 
affected (e.g. iTatulli et aT]|2007al) . To recalibrate the 
absolute continuum visibility level of the HR-K obser- 
vations, we make use of our earlier LR-HK observa- 
tions (Paper I) and the best-fit continuum model (Model 
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Fig. 2. — Spectra (top), visibilities (middle), and DPs (bottom) derived from our VLTI/AMBER observations with spectral resolution 
R = 12,000 obtained on 2010-03-03 (left) and 2010-08-26 (right). For the DP measurements, we add an arbitrary offset (dashed line) for 
clarity. 



DISK), which will be described in more detail in Sect. [31 
We investigated the influence of frame selection in the 
course of our data reduction procedure and find that the 
wavelength-differential signatures are very robust when 
selecting, for instance, the 80%, 50%, 20%, or 10% of 
frames with the best SNR. Given that the statistical 
noise increases with decreasing frame number and that 
low-SNR frames are down-weighted during the averaging 
process in a natural fashion, we decided to employ the 
observables from the full, unsclected data sets for our 
analysis. 

In order to associate the closure phase sign with the 
on-sky orientation, we us e a reference data s efl on the 
binary star 6 1 Orionis C (jKraus et al.1 l2009bl ) . To cal- 
ibrate the wavelength-differential phases, we make use 
of the fact that the continuum photocenter is systemati- 
cally offset with respect to the B^-line emission due to 
the presence of a close companion star. Using the well- 
established continuum closure phase sign calibration, we 
first determine the position angl(£3 (PA) of the compan- 
ion star (Sect. |3J and then adjust the DP sign in our 

2 The reference data set can be acc essed on the website 
http: / /www. stcfan-kraus.com/files/amber. htm 

J In this paper, all position angles are measured east of north. 



spectro-interfcrometric and spectro-astromctric data in 
order to match the direction of the continuum photo- 
center displacement. This procedure calibrates the DP 
sign unambiguously, and might serve also as a reference 
for other studies using AMBER's HR-mode. Recently, 
we employed this calibration in order to determine the 
rotatio n sense of the disk around the classical Be star 
/3CMi (jKraus et al.ll2012b[ ). Since /3CMi was observed 
using both AMBER's MR and HR-mode, we can also ex- 
tend the calibration to AMBER MR-mode observations, 
which resulted in a recali bration of the C Tau rotation 
sense (jKraus et al.1 l2012b| ) compared to earlier studies 
(IStefl et al.H2009h . 

The statistical error bars on the DPs vary strongly for 
the different baselines, which is both a result of the low 
fringe contrast at the longest baselines and of differences 
in the vibration properties of the UTs. These differences 
might give ovcrproportional weight to the short baselines 
in our x 2 -fitting proccsurc. Therefore, based on the typ- 
ical root-mean-square noise in the continuum channels, 
we include a minimum DP error of 20° for all baselines 
in the fitting procedure. 

2.2. VLT/CRIRES spectro-astrometry 
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Fig. 3.— VLT/CRIRES (i?. = 100,000) spectra obtained on 
V921 Sco in the Br7 and Mg II line. 



Spectro-astrometry uses high signal-to-noise (SNR) 
long-slit spectra to measure the centroid position of an 
unresolved object as function of wavelength. Since the 
centroid position can be measured with much hig her pre- 
cision than the size of the point-spread function (jBailevI 
11998b! ) , this method allows one to measure sub-mas pho- 
toccnter displacements in spectrally resolved emission 
lines. While the earliest astronomical applications in 
the field of star forma tion focused o n the detection of 
close companions (e.g. lBaiievlll998at ). later studies also 
successfully applied spectro-astrom etry to more general 
cases, such as outf low signatures (|Takami et al.l [20031 : 
iWhelan et al.l 120051 ) or the characte rization of the gas 
kinematics in protopla netary disks (Pontoppida n et al.l 
[20111: IGoto eFa!ll2012h . 

Our spectro-astrometric observations on V921 Sco were 
obtained using the VL T near-infrared hi gh-resolution 
spectrograph CRIRES (|Kaeufl et al.ll2004f ). The mea- 
surements were obtained on 2010-03-04, UT 07:51 to 
08:58 using an integration time of 5 s and a slit width of 
0.2", resulting in a spectral resolution of R = 100, 000. 
In order to apply the spectro-astrometry technique, we 
recorded spectra towards three different position angles 
(55°, 115°, 175°), and the corresponding anti-parallel po- 
sition angles (235°, 295°, 355°), complementing the po- 
sition angles probed by interfcrometry (Fig.Q]). To max- 
imize the flux in the slit, the spectra were recorded using 
the STRAP adaptive optics system, resulting in a typi- 
cal PSF Gaussian FWHM of ~ 115 mas. The measured 



spectra were corrected to the hclioccntric-baryccntric 
system using a heliocentric velocity of +28.36 kms -1 . 
The spectra were extracted using the ESO CRIRES data 
reduction pipeline (version 1.12.0). 

Besides the Br7 2.166078 /im line, which was covered 
by CRIRES detector #3, our observation also covered 
the Mg II 2.137 /im line in detecto r #1. This line has also 
been detected in other B[e] stars ([Clark et al.lll999[ ). To 
compute velocities in the local standard of rest (LSR), 
we assumed a systemic velocity of vq = —20 kms -1 , 
w hich we adopt from the f orbid den line measurements 
bv lBorges Fernandes et al.l (|2007l) . 

In order to correct the spectrum for telluric spectral 
features, we observed the early G-type star Hip84425 
and we modeled the intrinsic Br7-absorption using 
an infrared sola r spectrum recorded with ACE-FTS 
(|Hase et al.ll2010D . 

The final spectra are shown in Fig. [3] and reveal a 
very narrow (FWHM 0.0006560 /im or w 91 kms" 1 ) and 
double-peaked line profile with a very small peak separa- 
tion of just ~ 0.000184 fim (~ 25 kms -1 ). In addition to 
the Bi'7-line spectrum, which will be discussed and mod- 
eled in Sect.[H we have also recorded the Mg n 2.137 /im 
line. For the Mg n-line, the peak separation is signifi- 
cantly wider (~ 0.000500 /im or 70 kms -1 ), indicating 
that this line originates at smaller stellocentric radii in 
the circumstellar environment. This finding is consistent 
with our non-detection of a spectro-astrometric signal 
within this line. 

To derive the astrometric signal, we compute for 
each spectral channel y the beam centroid position 
([Pontoppidan et aI1l2011h in spatial direction 

X {V) - K EiFiiv) ' W 

where i defines the size of the virtual aperture (we use 
±1 pixel around the center of the PSF) and K is defined 
as the ratio between the flux in the virtual aperture and 
the total PSF (~ 1.5). The astrometric signals are de- 
rived for each position angle (X p ) and its anti-parallel 
counterpart (X a ) a nd then subtracted in o rder to remove 
potential artefacts (Branni gan et al.ll2006l ): 

X(v) = (X*(v)-X a {v))/2. (2) 

The interpretation of spectro-astrometric signals is 
complicated by the fact that this technique measures only 
the first-order momentum in the brightness distribution 
(relative position), but is insensitive to higher -order mo- 
ment a (angular size, asymmetry, kurtosis; see lLachaumd 
120031 ) as well as the continuum-emitting geometry, which 
limits this technique to relatively simple cases or requires 
additional model assumptions. In this study, we make 
use of the fact that photocenter displacements measured 
with spectro-astrometry are mathematically equivalent 
to the DPs measured in spectro-interferometry, which 
allows us to directly combine the spectro-intcrferometric 
(Visibility, DP, closure phases) and spectro-astrometric 
observables (photocenter displacements) for quantitative 
modeling fScct. l4~2")l . providing unique constraints on the 
AU-scalc spatial distribution and kinematics of the cir- 
cumstellar gas and dust. For this purpose, we translate 
the measured astrometric signal into the equivalent DP 
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using the relation 



2irX(v) 



(3) 



where a is the FWHM of the PSF measured in the spec- 
trum (typically ~ 150 mas). In order to give a similar 
weight to the AMBER and CRIRES measurements, we 
assume a minimum DP error of 1° for our CRIRES DPs. 

2.3. Magellan/FIRE near-infared spectroscopy 

We recorded a high-resolution, high-SNR near-infrared 
spectrum of V921 Sco using the FIRE Echelle spec- 
trograph (jSimcoe et al.l |2008( ) mounted at the Magcl- 
lan/Baade 6.5 m telescope. With a 0.45" slit, this in- 
strument provides a spectral resolution R = 8000 and 
a wide wavelength coverage from 0.8 to 2.5 /im. The 
spectrum was recorded on 2011-03-12 using an A-B-B- 
A dithering pattern with DITs of 10 s (high gain mode, 
to obtain high SNR in the z— , and J-band) and 5 s 
(low gain mode, in order to avoid saturation in the H- 
and AT-band). The slit was oriented along the disk po- 
lar axis derived from our interferometric observations 
(PA=56°). To correct for telluric features, we observed 
the AOV-type standard star HD 122945. The data re- 
duction was performed using the standard FIRE data 
reduction pipeline developed at MIT, which performs 
also a spectro-photometric calibration based using the 
photometry of the standard star. Trying different com- 
binations of target star / standard star data sets, we 
find that the derived spectral slopes are consistent on 
the level of a few percent, while the derived absolute flux 
levels show some larger scatter in the ~ 5 — 10% range. In 
order to improve on the calibration, we fit a fifth-order 
polynomical in order to extract the spectral slope and 
then use archival ISO photometry in order to recalibrate 
the absolute flux at 2.4/im which yields also a satisfy- 
ing match with 2MASS J-, H-, and A"-band photome- 
try (the derived flux densities match within 20%). The 
derived absolute-calibrated FIR E spectrum is corrected 
for reddening (Ay = 4.8 ± 0.2, iBorges Fernandes et al.1 
12001 using the extinction law from lCardelli et al.1 (jl989i 
R v = 3.1). 

At the longest wavelengths, the SNR in our final spec- 
trum is reduced compared to the J- and H-band since 
some of the recorded frames had to be rejected due to 
saturation in this part of the spectrum. In particular, 
the saturation occurs in sky lines, resulting in some nar- 
row spikes in the final corrected spectrum. 

3. RESULTS: CONTINUUM GEOMETRY 

In the following, we investigate the continuum struc- 
ture of the circumprimary disks using the VLTI/ AMBER 
low spectral dispersion data presented in Paper I. In this 
earlier study, the circumprimary disk emission was ap- 
proximated with a Gaussian brightness distribution (in 
the following denoted "GAUSS" model) , which provided 
a sufficient representation to extract the astrometric in- 
formation for the companion star. However, the rela- 
tively large Xr~ veL ^ ue °f 4-30 already indicates that this 
simple geometry does not provide an appropriate repre- 
sentation of the circumprimary disk structure. There- 
fore, we consider here a more realistic parameterization 
for the structure of the circumprimary disk, using a disk 



temperature-gradient model (in the following denoted 
"DISK" model). 

In the DISK model, the emission extends from an inner 
truncation radius i?; n to an outer radius i? ut and radi- 
ates as a black-body with T(r) = T- m (r / R ln )~ q , where 
Tj n denotes the temperature at the inner disk radius and 
q is the temperature power-law index. In addition, we 
introduce the inclination angle i and position angle tp, 
which define the projection ratio (cos i = i? m inor/-Rmajor) 
and on-sky orientation of the disk major axis. The in- 
clination is measured from the polar axis (i.e. i = 0° 
is pole-on). As in our earlier modeling attempts, the 
DISK model includes the photospheric emission of the 
primary and secondary star, which are parameterized 
by the companion separation (p) and PA (O), the an- 
gular extension of the circumsecondary material (given 
by a Gaussian with FWHM 8 b), and two parameters 
(■j=fk-(2 £tm), s) to describe the photospheric flux contri- 
butions of V921 Sco B to the total flux as function of 
wavelength: -j=^-(A) = -^-(2 /im) + s ■ (A — 2 fxm). 

We employ a Levenberg-Marquardt least square fitting 
procedure and find the best-fit solution by minimizing 
the likelihood-estimator \r — Xr.v+Xr <s>, where Xr.v an d 
Xr $ are the reduced least square between the measured 
and model visibilities and closure phases, respectively. 
The parameter uncertainties have been estimated using 
the bootstrapping technique. 

Given that the absolute calibration of the AM- 
BER+FINITO HR observations is not very reliable due 
to the use of long integration times and a potential 
residual phase jitter, we decided to recalibrate the abso- 
lute visibility level using the detailed continuum model 
discussed in Sect. [3] (Model DISK). Since the HR-data 
was recorded 2010, we extrapolated the orbital motion 
of the companion to this later epoch (p = 25.41 mas, 
9 = 339°), but note that the continuum geometry has 
only a marginal effect on our model-fitting results in the 
strong Br7-linc of V921 Sco. 

The DISK model provides a significant improvement 
compared to the GAUSS model, which also reflects in 
an improved ^-values of 1.86 (versus 4.88) for the 2008 
data and 3.02 (versus 3.39) for the 2009 data, respectively. 
In contrast to the GAUSS model, the DISK model also 
provides a physically motivated parameterization for the 
wavelength-dependent changes in the source geometry, 
allowing us to fit all wavelength channels simultaneously. 
The resulting best-fit parameters for the various models 
are listed in Tab. [TJ We plot in Fig. @] the measured 
observables versus the prediction of the DISK model. 
Besides the spatially resolved emission around the pri- 
mary component, we also checked whether the fit results 
can be improved by including spatially extended emis- 
sion around the northern (secondary) component, but 
were not able to find a significant improvement. 

4. RESULTS: Br 7 -LINE GEOMETRY 

4.1. 2-D photocenter analysis 

For a first, qualitative interpretation (Sect. I4.1j) of the 
gas kinematics, we derive the photocenter displacement 
of the line-emitting region with respect to the continuum 
emission using our CRIRES spectro-astrometric data. 
For this purpose, it is necessary to separate the line 
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Fig. 4. — VLTI/AMBER visibilities (LR-HK mode) of V921 Sco in the H- and if -band, plotted as function of spatial frequency. The 
different panels represent position angle bins, covering 10° each. The grey model curves correspond to our best-fit temperature gradient 
disk model (DISK) and were computed for the lower (1.45 pm, dotted line), central (2 pm, dashed line), and upper (2.55 pm, solid line) 
part of the wavelength range. For clarity, the data points are color-coded according to the wavelength of the corresponding spectral channel 
(blue: 1.4 < A < 1.9 /im; green: 1.9 < A < 2.15 pm; red: 2.15 < A < 2.5 fim). 



TABLE 1 

Model-fitting results for the VLTI/AMBER continuum observations (Sect. 01 



Secondary star Circumprimary disk 

Epoch Model Q p 8 B s 8 A i <P Rin Rout q Xr,v X?,* X? 

[°] [mas] [mas] (at 2 pro) [jim -1 ] [mas] [°] [°] [mas] [mas] 



2008 


GAUSS 


353.8 


25.0 


< 0.2 


0.054 


-0.0056 


7.5 


50.3 


147.8 








4.30 


3.09 


4.88 






±1.6 


±0.8 




±0.018 


±0.019 


±0.2 


±1.9 


±4.3 














2008 


DISK 


354.0 


24.9 


< 0.7 


0.066 


-0.051 




48.8 


145.0 


1.59 


> 9.7 


0.36 


1.86 


1.86 


1.86 






±1.0 


±0.5 




±0.04 


±0.01 




±4 


±8.4 


±0.25 




±0.08 








2009 


GAUSS 


347.3 


25.5 


< 0.2 a 


0.054" 


-0.0056" 


7.5 a 


50. 3 a 


147.8" 








4.17 


1.78 


3.39 






±1.0 


±1.2 


























2009 


DISK 


347.4 


25.4 


0.7 Q 


0.086 


-0.051 Q 




48.8 a 


145.0° 


1.59 Q 


9.7 a 


0.36 a 


3.26 


2.52 


3.02 






±1.1 


±0.7 




±0.012 























a In our fitting procedure, this parameter was kept fixed. 



spectro-astrometric signal from the underlying contin- 
uum contributions. To determine the continuum level, 
we fit a high-order polynomial function X c to the astro- 
metric signal derived from the continuum channels. This 
function is then subtracted from the astrometric signal in 
the line channels and weighted by the continuum-to-linc 
flux ratio F c /Ff. 

= (x(v) - x c(v)) (1 + fc(w)/*K«)) (4) 
The derived photoccnter of the line velocity channels 



is significantly offset with respect to each other (clearly 
indicating the gas kinematics) , and also shows a displace- 
ment with respect to the photocenter of the continuum 
channels (Fig. [6]). This displacement reflects the fact that 
the center of gravity of the continuum emission does not 
coincide with the location of the primary star, but is 
displaced towards the position of the companion. We 
calculate the center of gravity both for the line channels 
and the continuum channels and compute their relative 
displacement vector v c -i, which is related to the com- 
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Fig. 5. — VLTI/AMBER closure phases and model closure phases for our best-fit model DISK. In each panel, we give the projected 
length and PA for two of the three employed baselines (the values for the third baseline are given by the closure relation). 



panion separation p, position angle 0, and the flux ratio 
#^by 

■rtot J 

F B ( sine \ , . 

= ^ Use )■ (5) 

Based on our continuum model fits from 2008 and 2009 
(Tab. [1]), we fix p = 25.13 mas and then determine the 
companion position angle 6 and the flux ratio Fb /-Ftot 
from the CRIRES observations. Applying this procedure 
yields & = 353 ± 9° (which is in reasonable agreement 
with the astrometry determined with AMBER, Tab. [lj 
and Fs/Ftot = 0.053, which is lower than the flux ra- 
tio determined with interferometry for epochs 2008 and 
2009, possibly indicating variability. 

The photocenters in the blue- and red-shifted line 
wings arc displaced in opposite direction with respect 
to each other. Also, the highest gas velocities (dark 
red points and dark blue points in Fig. [S]) emerge from 
smaller stellocentric distances than intermediate veloci- 
ties, which suggests a rotation-dominated velocity pro- 
file. From the distribution of the individual photo- 
center offsets, we derive the disk plane orientation to 
153.2 ±0.6°. 

In contrast to CRIRES, AMBER is able to spatially re- 
solve the geometry of the line-emitting region, entering a 
regime where higher-order geometric effects are probed. 
Therefore, in the following two sections, we will employ a 
quantitative modeling in order to interpret the combined 
spectro-astrometric and spectro-intcrfcrometric observa- 
tions. 



TABLE 2 

Model-fitting results for the VLTI/AMBER (HR-K) and 
VLT/CRIRES Br7-LINE data (Sect. [p. 



Parameter 






Value 


Inner emission radius 


Rin 


[mas] 


$ 0.5 


Outer emission radius 


Rout 


[mas] 


5.8 ±0.3 


Position angle 





H 


145.0 ±1.3 


Inclination 


i 


H 


54 ±2 


Stellar mass 




[M ] 


(5.4 ±0.4) ■ (d/1150 pc) 


Radial intensity index 


q 




-1.0 ±0.1 


Opacity index 


K 




0.5 ±0.2 



4.2. Model fitting: Keplerian disk 

Given the indications for a rotation-dominated veloc- 
ity field provided by our model-independent photocenter 
analysis (Sect. l4~Tj) . we first test whether a Keplerian ve- 
locity field might reproduce our spectroscopic, spectro- 
astrometric and spectro-interferometric data quantita- 
tively. 

For this purpose, we employ a Keplerian disk model 
(v(r) = \/GM+/r), which we have already successfully 
applied to AMBER high sp ectral dispersion d ata on 
the classical Be star /3CMi (jKraus et all l2012b[ ). The 
model assumes that the line-emitting gas is optically thin 
and located in a geometrically thin disk, where the ra- 
dial intensity profile is parameterized with a power-law 
(Ii{r) oc r q ). The gas emission extends from an in- 
ner truncation radius (R- m ) to an outer truncation ra- 
dius (i?out) ; which we parameterize with a Fermi- type 
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velocity LSR [km/s] 



-150 -100 -50 50 100 150 




E 
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115° 



175° 



CO 

o 

LU 
Q 

"D 




0.5 -0.5 -1 
dRA [mas] 



2.164 2.165 2.166 2.167 

Wavelength [fim] 

Fig. 6. — Left: Spectrum (top) and spectro-astrometric signal (2nd-4th panel from top) derived from our VLT/CRIRES (R = 100, 000) 
observations. Right: Derived 2-D photocenter displacements, where the different spectral channels are represented by different colors, 
matching the color coding in the left panel. From the photocenter displacements measured in the different spectral channels, we derive the 
disk rotation plane to 153.2 ± 0.6°. The continuum spectral channels (grey points) arc offset with respect to the center of light of the line 
emission (origin of the coordinate system) due to the presence of the companion. 



function in order t o avoid artifical edges (see Eq. 1 in 
IKraus et al.1 12008al where the width of the truncation 
region e was chose n to 0.1). In our study on j3 CMi 
(|Kraus et al.ll2012bfl . we investigated the DP signatures 
induced by photospheric absorption and found that even 
for a relatively extended stellar surface (equatorial radius 
of 0.36 mas), the induced DPs are < 0.5°. For V921 Sco, 
the influence of photospheric absorption is likely ~ 50- 
times smaller due to the smaller apparent stellar ra- 
dius (~ 0.07 mas) and the significantly larger equivalent 
width of the Bi'7 emission line. Therefore, these signa- 
tures are about four orders of magnitude smaller than the 
disk kincmatical signatures and can be safely neglected 
in our modeling process. 

As discussed in Sect. 14.11 the elliptical distribution of 
the derived photocenter vectors might indicate opacity 
effects, which cause the more distant parts of the disk to 
appear fainter than the disk parts facing the observer. 



In order to include this effect in our model, we assume 
that the disk is embedded in an medium of constant den- 
sity. After a distance x traveled, the emitted intensity 
Iq is then reduced to I(x) = io ex p(~ K2 0- Combin- 
ing this line-emission model with our best-fit continuum 
disk+companion model (DISK, Tab.Q]) allows us to pro- 
duce model channel maps (Fig. bottom) , from which 
we compute line profiles, visibilities, and DPs for com- 
parison with our data (Fig. top). 

As line broadening mechanisms, we include therm al 
and turbulent Dopplcr broadening (jPietu et alj |2007) 



AV 



2kT -, 



J turb' 



(6) 



where fturb is the turbulent velocity, which is typically 
negligible for He rbig Ae/Be stars (< 0.5 kms -1 , e.g. 
iPietu et al.ll2007l ). In the disk surface layer, the gas tem- 
perature T gas will be significantly higher than the dust 



V921 Scorpii: Geometry and kinematics of the circumprimary disk on sub-AU scales 



9 



velocity [km/s] velocity [km/s] velocity [km/s] 



-300-200-100 100 200 300 -300-200-100 100 200 300 -150 -100 -50 50 100 150 




2.162 2.164 2.166 2.168 2.17 2.162 2.164 2.166 2.168 2.17 2.164 2.166 

Wavelength Qim] Wavelength [|im] Wavelength [|im] 




2.5 0.0 -2.5 2.5 0.0 -2.5 2.5 0.0 -2.5 2.5 0.0 -2.5 2.5 0.0 -2.5 

dRA [mas] dRA [mas] dRA [mas] dRA [mas] dRA [mas] 



Fig. 7. — Top: Spectra (1st row), visibilities (2nd row), an d dif ferential phases (3rd row) derived from our AMBER and CRIRES Br7 
observations, overplotted with the Keplerian disk model ('Sect, |4T2)| , For each interferometric baseline, wc give the corresponding projected 
baseline length and PA in the label. In order to make the visibility profiles in the left panel clearly distinguishable, we added an offset of 
0.2 to the visibility corresponding to the baseline 46.6m/22°. Bottom: To illustrate our kinematical model we show channel maps for some 
representative wavelengths, including the Br7-linc and the continuum emission. The continuum-emitting disk is truncated at a distance 
of 1.59 mas, while the Br7-emitting material extends inwards to a few stellar radius (< 0.5 mas). The companion star V921 Sco B is also 
included in these simulations, but located outside of the FOV in these model images. 



temperature at a given stellocentric radius. We assume 
T gas = 5 Td U st, as suggested by the therma lly decou- 
pled r adiative-hydrodynamics simulations by iThi et aLl 
(|201lD . but note that the precise value affects mainly the 
line width and does not affect our general conclusions 
(see also Sect. 15. f I) . 

Free parameters in our model are the mass of the cen- 
tral star Ma, the inner and outer disk radius (R nl , R ut), 
the radial intensity power-law index q, the disk inclina- 
tion i, and the opacity index k. 

We vary these parameters systematically on a param- 
eter grid and select the best-fit model with the best 

Xr = Xly + Xr,<j> + Xr,Fl where \r,4, and Xr,F are the 

reduced chi-squared between the model and measured 
different phase and spectrum, respectively. 

Remarkably, this simple disk model can reproduce im- 
portant features in our data, in particular: 

(a) Our model can reproduce the double-peaked By- 
line profile (Fig. [7J 1st panel) reasonably well. One 
line profile characteristic, which is not reproduced 
by our model concerns the weak asymmetry ob- 
served in Br7-line profile, with a slightly stronger 



red-shifted line wing. 

(b) Our AMBER observations (Fig. [JJ 2nd panel) re- 
veal M-shaped visibility profiles on two of three 
baselines (36.7m/61° and 53.4m/89°), indicating 
that the angular extension of the line-emitting re- 
gion increases for low gas velocities. This effect 
is expected for a rotation-dominated velocity field, 
where the azimuthal velocity decreases as function 
of radius and the visibility profile shape is reason- 
ably well reproduced by our model for these base- 
lines. On the third baseline (46.6m/22°), no cen- 
tral visibility drop has been observed. We suspect 
that this effect is related to more subtle radiative 
transfer or line broadening effects, which are not 
included in our simplistic model. The visibilities 
at low gas velocities (i.e. in the line center) would 
be most sensitive to such small-order effects and wc 
leave it to future studies using full radiative trans- 
fer modeling to investigate this in more detail. 

(c) The measured AMBER and CRIRES DPs (Fig. [7J 
3rd panel) show very interesting signatures, includ- 
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ing asymmetric S-shapcd signatures and V-shaped 
signatures, which appear sometimes in the bluc- 
and sometimes in the red-shifted line wing. These 
signatures constrain the disk kinematics primarily, 
but contain also contributions from the continuum 
photocentcr displacement caused by the compan- 
ion star. Our combined binary star plus Keplerian 
disk rotation model reproduces the large variety of 
signatures reasonably well. The strongest residu- 
als between the model and the data are observed 
at the longest AMBER baselines, where the data 
is extremely sensitive to small-order kinematical ef- 
fects, but where the SNR is also reduced due to the 
lower fringe contrast. 

The parameters corresponding to the best-fit model are 
listed in Tab. [2] and will be discussed in more detail in 
Sect. El 

We also investigated whether a significant fraction of 
the Br7 emission might be associated with the secondary 
star (V921 Sco B) instead of the circumprimary disk, pos- 
sibly indicating active accretion from the circumbinary 
disk onto the secondary star. Considering the strong 
measured displacement of the line photocenter in the di- 
rection towards the primary (Fig. H3 right) and the good 
quantitative agreement of the displacement vector with 
the flux- weighted binary astrometry vector (Sect. 14. ip . 
it is clear that the secondary is not the dominant Br7- 
emitting component. In order to further quantify these 
constraints, we introduced the fraction of Br7 emission 
associated with the secondary star to the total Br7-flux 
(Fq / Fto^Br-y as an additional parameter in our kinemat- 
ical model and find that (FB/F to t)Br7 < 0.11. 

5. INTERPRETATION 
5.1. Spectroscopy 

Our FIRE observations (Fig. |5|) reveal a very rich 
near-infrared spectrum. Simil ar to the optical spectrum 
(|Borges Fernandes et al.1l2007f ). the near-infrared regime 
is dominated by strong hydrogen line emission. As listed 
in Tab. 02 we detect 27 lines from the Brackett series 
(Br7 to Br29), 24 lines from the Paschen series (Pact to 
Pa24), and at least 10 lines from the Pfund scries (Pfl4- 
24) . Given that the recombination physics of hydrogen 
is well known, one can use these line decrements in order 
to derive information about physical conditions of the 
line-emitting gas. 

In order to model the line flux in the optically thin 
appro ximation, we extracted from iStorev fc Hummer! 
(|1995| ) the line emissivities e u i for Case B recombination 
for all Balmer (I ~ 3), Brackett (I = 4), and Paschen 
(I = 5) hydrogen line transitions (u — I) with u < 25. 
The line emissivities include collisional transitions, which 
makes our computation applicable even for relatively 
high gas densities. Assuming that the emitting gas is 
isothermal (T = 10, 000 K) and uniformly distributed in 
the volume V, the received line flux at distance d is then 
given by 



where N e and N p are the electron and proton den- 
sity per unit volume. In order to approximate the 



emitting volume V, we estimate from our Br7 spectro- 
intcrfcromctric data the characteristic stellocentric emis- 
sion radius. For this purpose, we generate a velocity- 
integrated line map and measure the half-light emission 
radius to ?'i; nc = 3.0 mas=3.5 AU (for d = 1.15 kpc). 
Assuming that the line-emitting gas is located in a hot 
disk surface layer with a characteristic scale height h of 
h/ry mc = 0.2, as suggested by th e thermally decoupled 
simulations from iThi et al.l (|2011[ ). the emitting volume 
is given by V = nr^^ ■ (h/r\- ulc ). We employ an iter- 
ative procedure, where we start by assuming a low gas 
density, which is then adjusted in order to reproduce the 
measured line fluxes. 

Then we recompute the line emissivities and repeated 
the computation until convergence is reached. As shown 
in Fig. [51 we require rather high gas densities between 
N e = N p = 2 x 10 19 m~ 3 and 6 x 10 19 m" 3 , which sup- 
ports the scenario that the gas is located in a dense disk 
instead of a low-density halo. For the high-level transi- 
tions (u > 16, i.e. Pal4-24, Brl3-29, and Pfl2-23), the 
line fluxes can be reproduced well with a single gas den- 
sity (N e = 6 x 10 19 m~ 3 ), suggesting that the optical 
thin approximation is well justified. The low-level tran- 
sitions (u <= 16, i.e. Pa/3-Pal3, Br7-Brl2) and more 
consistent with densities around N e = 2 x 10 19 m~ 3 , 
which might indicating that these transitions originate 
from more extended disk regions, where the disk surface 
density is lower. We would like to note that the derived 
densities depend to some extend also on the assumed 
gas temperature. Varying the temperature from the as- 
sumed T — 10, 000 K to other realistic values (e.g. 5,000- 
20,000 K), would change the derived densities by a factor 
of ~ 2 — 3. Based on our simplistic modeling, the derived 
density should be treated as an ordcr-of-magnitude esti- 
mate and future radiation-hydrodynamic modeling will 
be required in order to fully exploit the rich informa- 
tion provided by our combined spectro-interferometric 
and spectroscopic dataset. 

It is interesting to compare the derived electron den- 
sity to the values predicted by state-of-the-art disk mod- 
els. For instance, for a typical disk around a Herbig Ae 
star (e.g. AB Aurigae, 2.4 Mq) with a gas surface den- 
sity of 10 4 kgm~ 2 , the predicted midplane gas density 
is ~ 6 x 10 17 m" 3 (jDullemond fc Monnier1l2010f ) at the 
dust sublimation radius (1.1 AU). Our density measure- 
ment of N e = 2. ..6 x 10 19 m~ 3 probes similar spatial 
scales in the V921 Sco disk and is one to two orders of 
magnitude higher than this value, suggesting that the 
disk around V921Sco is exceptio nally massive. This re - 
sult is in line with the study by iHenning et all (|1998f) . 
who determined the total gas mass in the millimeter 
cores around 25 Herbig Ac/Be and FU Orionis stars 
and measured the highest total gas mass in the core 
around V921 Sco (M gas = 40 M Q ; core size: 27 x 27" or 
~ 30,000 x 30,000 AU at 1.15 kpc). Observations with 
our spectroscopic & spectro-interferometric approach on 
a larger sample of protoplanetary disks will be necessary 
in order determine whether the derived high gas density 
indicates an intrinsic property of the V921 Sco disk (pos- 
sibly due to the young age of the source) or simply reflects 
the predicted disk surface density scaling law with stellar 
mass (S oc M cx M 2 . ICalvet et~all[200l . 

Besides hydrogen recombination lines, other identified 
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lines include Fe II, [Fe II], C I. He I, O I, N I. Mg II, 
Al II, w hich have also be en identified in the B[el stars 
CI C am (IClark et al.lll999h and H P 50138 (iJaschek et al l 
1991. or in the LBV rj Car (jDamineli et all 119981 : 
Smith fe Davidson|[200TI) . Besides the total of 90 identi- 



TABLE 3 

Lines identified in our V921 Sco FIRE 

SPECTRUM AND MEASURED LINE EQUIVALENT 
WIDTHS (EW) AND INTEGRATED LINE FLUXES F 
(FOR LINES IN THE J-BAND). 



ficd lines, we detect in our spectra a few dozen additional, 
yet unidentified lines and we encourage detailed follow- 
up spectroscopic studies in order to investigate the rich 
chemistry of this source. 

Our wavelength range also covers the CO first overtone 
bandheads between 2.3 and 2.4 /urn, which allows us to 
search for the CO first overtone bandheads, which we de- 
tect neither in absorp tion nor emis sion. This is particu- 
larly interesting, since lKra"ui (|2009l ) suggested that these 
lines provide a good diagnostic tool to distinguish be- 
tween a pre- and post-main-sequence evolutionary phase 
in B[e] stars. According to her computation, observable 
13 CO bandhcad emission can only be produced in evolved 
stars, but should be absent in pre- main-sequence B[e] 
stars. 

Finally, we investigated whether we find indications 
for variability of the Br7-line during the ~ 6 months cov- 
ered by our CRIRES and AMBER observations. For this 
purpose, we convolved the CRIRES and AMBER HR-K 
data to the same spectral resolution and overplot the line 
profiles in Fig. I10| finding no indications for variability, 
neither in line strength nor line profile. 

5.2. The circumprimary disk 

Using a temperature-gradient disk model (DISK), 
we determine the inner continuum disk radius i?; n to 
1.59 mas. We can compare this value with the expected 
location of the dust sublimation radius, using the lumi- 
nosity and distance estimate of iBorges Fernandes et al.l 
(|2007h . L+ = (10 ± 3) x 10 3 L Q and d = 1150 ± 150 pc. 
Assuming grey dust opacities and a standard dust sub- 
limation temperature of T su bi = 1500 K, the expected 
dust sublimation radius (including the effect of back- 

warming from the dis k) is i? su bi = y L± /4naT^ nhl 
(jDullemond et all I2001D , where a denotes the Stcfan- 
Boltzmann constant. With the given large luminos- 
ity and distance uncertainties, this corresponds to a 
value of 6.9 ± 1.1 AU or 6.0 ± 1.9 mas, which is con- 
siderably larger than the measured inner disk radius 
of i?i n = 1-59 mas. Assuming higher dust sublima- 
tion temperatures of 2000 K (i? su bi = 3.9 ± 0.6 AU or 
3.4±1.0 mas) also does not solve this discrepancy. There- 
fore, we conclude that V921 Sco belongs to the group 
of Herbig Be disks that are " undersized" with respect 
to th e size-luminosity relation dMonnier fe Millan-Gabet] 
12002). as already discu ssed bv iKraus et all (|2008bf ) 
and iKreplin et all ()2012l ). Different physical scenar- 
ios have been proposed in order to explain this ef- 
fect, including gas absorption ([Monnier fe Millan-Gabetl 
[200a . the emission f rom gas located in s ide of the dust 
sublimation radius (lEisner et al.l 120041: iMonnier et al.l 
[2005t IKraus et al.l l2008at iTannirkulam et al .1 120081 ). or 
the presence of a h ighly refractive dust grain species 
(jBenistv et al.ll2010( ). Multi- wavelength interferometric 
studies with even longer baseline lengths as well as de- 
tailed gas and dust radiative transfer modeling will be 
required in order to ultimately settle this question. 
We determine the outer disk radius to ~ 9.7 mas. How- 



Wavelength 


Line 


EW 


F 


H 




[nm] 


[10- 16 Wm- 2 ] 


0.8304 


Pa24 


-0.053 a 


14.7 


0.8312 


Pa23 


-0.073 


20.1 


0.8321 


Pa22 


-0.079 


21.6 


0.8332 


Pa21 


-0.087 


23.7 


0.8343 


Pa20 


-0.110 


29.7 


0.8357 


Pal9 


-0.153 


41.0 


0.8372 


Pal8 


-0.175 


46.4 


0.8390 


Pal 7 


-0.199 


52.2 


0.8411 


Pal6 


-0.242 


62.7 


0.8436 


Pal5 


-0.300 


76.6 


0.8447 


O I 


-1.343 


340.4 


0.8465 


Pal4 


-0.369 


92.6 


0.8500 


Pal3 


-0.443 


108.9 


0.8543 


Pal2 


-0.432 


103.7 


0.8596 


Pall 


-0.339 


79.1 


0.8663 


PalO 


-0.394 


88.8 


0.8748 


Pa9 


-0.348 


75.3 


0.8820 


O I 


-0.012 


2.5 


0.8860 


Pa8 


-0.447 


92.2 


0.8927 


Fe II 


-0.097 


19.5 


0.9013 


Pa7 


-0.462 


89.9 


0.9127 


Fe II 


-0.074 


13.9 


0.9177 


Fe II 


-0.169 


31.2 


0.9202 


Fe II 


-0.121 


22.1 


0.9227 


Pa6 


-0.500 a 


90.8 


0.9256 


Fe II 


-0.080 a 


14.4 


0.9261 


O I 


-0.070 a 


12.6 


0.9392 


M I 


-0.071 a 


12.3 


0.9406 


C I 


-0.060 


10.4 


0.9543 


Pae 


-0.660 


110.1 


0.9893 


Fe II 


-0.045 


6.9 


0.9910 


Fe II 


-0.029 


4.5 


0.9956 


Fe II 


-0.070 


10.7 


0.9997 


Fe II 


-0.624 


94.4 


1.0047 


Pa<5 


-1.132 


169.9 


1.0108 


Al II 


-0.043 a 


6.4 


1.0174 


Fe II 


-0.143 


21.0 


1.0501 


Fe 11 


-0.549 


76.6 


1.0688 


C I 


-0.163 


22.2 


1.0829 


He I 


-0.760 


101.7 


1.0862 


Fe II 


-0.503 


67.1 


1.0917 


He I 


-0.207 


27.4 


1.0935 


Pa7 


-1.793 a 


237.2 


1.1127 


Fe II 


-0.468 


60.6 


1.1290 


O I 


-1.874 a 


238.9 


1.1755 


C I 


-0.117 


14.3 


1.2467 


N I 


-0.083 a 


9.6 


1.2815 


Pa/3 


-4.030 a 


457.5 


1.3166 


O I 


-0.321 


35.6 



a These lines suffer from line-blending, which intro- 
duces some uncertainty to the derived EW values. 

Due to the high if -band continuum flux and line 
brightness, these lines likely reached the non-linearity 
regime of the FIRE detector, resulting in an underes- 
timation of the real EW. 

ever, this value should be treated as a lower limit, since 
colder parts of the disk might extend to larger stellocen- 
tric radii. We can compare this derived lower limit with 
theoretical pred ictions for the disk truncatio n radius in 
binary systems. lArtvmowicz fe Lubowl (|1994ft computed 
the resonance- induced truncation radius for circumstellar 
and circumbinary disks and predicts the circumprimary 
disk to be truncated at r/a <~ 40... 20% of the orbit ma- 
jor axis (for a binary mass ratio of 0.3), depending on 
the orbit eccentricity and disk viscosity parameter. This 
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Fig. 8. — Near-infrared J-band (fop), H-band (middle), and A'-band (bottom) spectrum recorded with the Magellan/FIRE spectrograph. 
The spectrum is dominated by strong hydrogen emission lines from the Paschen, Brackett, and Pfund series (grey lines) and metallic lines 
from Fe II (red, dashed lines), [Fe II] (red, dotted lines), Mg II, N I, Al II, C I, O I (green lines). A complete list of the identified lines can 
be found in Tab. [3] For the Pfund lines, we show besides the clearly identified transitions (Pfl4-24, grey dashed lines) also the higher-order 
line transitions (Pf25-100, grey dotted lines) in order to mark the location of the Pfund discontinuity. 



range agrees reasonably well with our derived value of 
Rout/p — 39%. For a detailed comparison, a full orbit 
determination and mid-infrared and sub-millimeter in- 
tcrfcromctric observations of the colder material in the 
outer regions of the circumprimary disk and in the cir- 

cumbinary gas and dust reservoir will be requi red. 

Both theore ti cal (e.g. iNatta et all 120011 : 
iDullemond eTaTl [20011 
observational studies (M o 
20092 iBenistv et all 120111 



llsella fc Nattal 
mier et alJl2006t 



20051) and 



Kraus et al.1 



suggest that the inner dust 
rim in protoplanetary disks has a complicated, possibly 
puffed-up vertical structure, which can result in an 
asymmetric brightness distribution. Intcrferomctric 
data, such as presented in this study, is able to detect 
these asymmetries and to constrain the corresponding 
disk models, in particular using phase-closure capabili- 
ties. However, it is clear that the closure phases in our 
AMBER data are dominated by the signatures of the 
newly-detected companion, while potential asymmetries 
due to the vertical disk structure would appear only 



as secondary effects. Such signatures would appear 
in particular at long baseline lengths and close to 
visibility minima, where we also observe some significant 
deviations from our simple flat disk+companion model 
(Fig. [5]). However, it would be difficult to constrain 
these more complex physical models with the current 
uncertainties on the stellar parameters and the lim- 
ited available uw-coverage. Therefore, we leave it to 
future investigations to better constrain the inner dust 
disk geometry around V921Sco A. Alternatively, the 
remaining residuals might indicate the orbital motion 
of the companion within the five months during which 
the 2008 AMBER LR-HK data has been recorded (sec 
Tab 1 in Paper I). 

Using our kincmatical modeling, we show that the ve- 
locity field in the disk is in Keplerian rotation (Sect. |4~2"|) 
and that the line-emission emerges from a relatively high- 
density region (N e — 2. ..6 x 10 19 m -3 ; Scct. l5.1| ). likely in 
a hot disk surface layer. Between the high-level (u > 16, 
i.e. Pal4-24, Brl3-29, and Pfl2-23) and low-level tran- 
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Fig. 10.— VLT/CRIRES and VLTI/AMBER Br7-linc spectra 
measured in March and August 2010, which have been convolved 
to a common spectral resolution of R = 12, 000. 



sitions (u < 16, i.e. Pa/?-Pal3, Br7-Brl2), we detect a 
change of slope in the Paschen/Brackett decrements, in- 
dicating that the high-level transitions emerge from high- 
density regions closer to the star than the low-level tran- 
sitions. In our best- fit model (Tab. [21 Fig.[7J), the orienta- 
tion of the opacity screen is such that the south-western 
part of the disk appear brighter, which suggests that the 
north-eastern disk axis is facing towards us (which is also 
consistent with the orientation of the large-scale bipolar 
nebula, Fig. [11] A-C). Arguably the most surprising pa- 



rameter value in our best-fit model concerns the derived 
stellar mass. For an assumed distance of 1.15 kpc, we 
yield best agreement with 5.4 ± 0.4 M Q , which is hardly 
consistent with the spect ral classification^! and the pre - 
diction mass of 9±1 M Q (jBorges Fernandas et al.|[2007h . 

In order to explain this surprising result, we propose 
two plausible scenarios: 



(1) 



The value of 1 .15 ± 0.15 kpc 
(jBorges Fernandes et al.l 120071) might under- 
estimate the real distance to V921 Sco. This 
distance has been determined based on the equiv- 
alent width of the Na I and Ca II-K interstellar 
absorption lines and therefore relies on a sta- 
tistically es t ablishe d calibration. In particular, 
lLopes et al.l (|1992t) employed the same method 
(also on the Na I line) and determined the distance 
to 2.5 kpc. A di stance estimate ^ 2 kp c has also 
been obtained by iMcGregor et a l. (|1988f l based on 
optical/near- infrared spectroscopy and photome- 
try. Adopting such a larger distance would result 
in a stellar mass of 12 M Q (for 2.5 kpc) or 9 Mq 
(for 2 kpc), which would be more consistent with 
the spectral classification. 



4 It should be noted that the earlier spectroscopic studies were 
not informed about the binarity of V921 Sco, which has likely also 
resulted in a somewhat biased spectral classifications of this in- 
triguing source. However, given the extremely rich recombination 
line spectrum and the presence of He I emission lines, it seems safe 
to assume that the ionizing central source should be of at least 
B-type. 
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TABLE 4 

Continued: Table [3] for lines in the H- and 

H -BAND . 



Wavelength 


Line 


EW 


F 


M 




[nm] 


[10- 16 Wm" 2 ] 


1 A "7CIQ 

1.4 /yo 




-U.Uol 


Q 1 

o.X 


1.4811 


Br28 


-U.UoU 


3.0 


1.4827 


Br27 


-U.U4Z 


4 . 2 


1.4844 


Br26 


-0.051 


5.1 


1.4863 


Br25 


-0.057 


5.7 


1.4884 


Br24 


-0.072 


7.3 


1.4907 


Br23 


-0.075 a 


7.5 


1.4934 


Br22 


-0.096 


9.6 


1.4963 


Br21 


-0.117 


11.7 


1.4997 


Br20 


-0.141 


14.1 


1.5035 


Brl9 


-0.161 


16.1 


l.OU ( o 


D r 1 O 


n iQi 
-u. lyi 


iy . i 


i.oi/y 


±31 1 ( 


n 01 o 

-U.Z1Z 


0~\ 1 
z 1 . 1 


l.Oloo 


r>rlo 


-U.ZOU 


OA Q 


l.OZOO 


r>rlO 


-U.z / y 


07 ft 


l.OOOO 


15114 


-u.oyo 


oo.y 


1.04D0 


J3llO 


-U.4oU 


AO 1 
4Z. 1 


1.000Z 


lorlz 


-U.0Z4 


01. U 


±.ooyo 


Rr1 1 


-U.044 


OZ.D 


1.0 ( ou 


EVi TT 

rC 11 


-U.1U0 


in 1 
1U. 1 


l.Oo (D 


131 1U 


-U.0 (U 


04. ( 


1.6105 


Bi9 


-0.631 


59.9 


1.6403 


Br8 


-0.667 a 


62.5 


1.6435 


[Fe III 

L J 


-0.019 a 


1.8 


1.6769 


[Fe II] 


-0.076 a 


7.0 


1.6787 


Fe II 


-0.046 a 


4.2 


1.6802 


Br7 


-0.719 a 


66.2 


1.6873 


Fe II 


-0.309 


28.4 


1.6894 


Fe II 


-0.015 a 


1.4 


1.7006 


He I 


-0.032 a 


2.9 


1.7111 


[Fe II] 


-0.013 


1.2 


1.7357 


Br6 


-0.521 


47.1 


1.7414 


Fe II 


-0.102 


9.2 


1.7451 


[Fe II] 


-0.054 a 


4.9 


1.9440 


Br<5 


-0.808 b 


70.3 


1.9738 


Fe II 


-0.019 


1.7 


2.0590 


Mg II 


-0.219 b 


19.1 




Mg II 


-U.UO-L 


i. ■ ) 


2.1447 


Mg II 


-0.018 


1.6 


2.1649 


Br7 


-0.705 b 


62.5 


2.3532 


Pf23 


-0.028 


2.5 


2.3591 


Pf22 


-0.036 


3.2 


2.3657 


Pf21 


-0.036 


3.2 


2.3731 


Pf20 


-0.054 


4.9 


2.3815 


Pfl9 


-0.064 


5.8 


2.3912 


Pfl8 


-0.078 


7.0 


2.4023 


Pfl7 


-0.075 


6.7 


2.4151 


Pfl6 


-0.092 


8.2 


2.4300 


Pfl5 


-0.114 


10.2 


2.4477 


Pfl4 


-0.162 


14.4 


2.4687 


Pfl3 


-0.245 


21.5 



a These lines suffer from line-blending, which intro- 
duces some uncertainty to the derived EW values. 

Due to the high i^-band continuum flux and line 
brightness, these lines likely reached the non-linearity 
regime of the FIRE detector, resulting in an underes- 
timation of the real EW. 



(2) The disk might not be associated with the massive 
(early B-type) star in the system, but orbit around 
the intermediate-mass (late B-type) companion in- 
stead. This would explain the relatively low de- 
rived mass, while the massive primary star could 
still provide sufficient ultraviolet flux to ionize the 
disk surface layer and the ambient low-density ma- 
terial resulting in the forbidden line-emission. In 
this scenario, the accretion would occur from the 



circumbinary disk onto the disk around the late B- 
type star. The massive primary would rest closer 
to the center-of-mass and be effectively shielded 
from mass-infall by the orbiting companion. A sim- 
ilar dynamical configurations seems also plausible 
for other Herbig B[e]-star binary systems such as 
MWC 361 A, where the active accretion is clearly 
associat ed with the less massive component in t he 
system (jMonnier et al.ll2006l: lAlecian et al.ll2008t ). 

In order to distinguish between these scenarios, it 
will be necessary to accurately measure the flux ratio 
of the two stars as function of wavelength. Our cur- 
rent measurement yield similar flux ratios between the 
two stars in the H- and JT-band, with some indications 
that V921Sco A (i.e. the star associated with the cir- 
cumstellar material) is of earlier spectral type (Paper I), 
favouring scenario (1). However, it should be noted that 
this color measurement is only significant at the 2a-level 
((F B /F A ) H = 0.83 ±0.15, (F B /F A ) K = 1.18 ±0.12, Pa- 
per I) and might be biased, for instance in case of in- 
homogeneous extinction towards the two stars. In both 
scenarios, it is clear that a careful re-evaluation of the 
spectral classification based on the spectral energy dis- 
tribution and the spectroscopic diagnostics will be nec- 
essary, which, is out of the scope of this paper. 

Independent of these uncertainties, we can derive new 
insights on the evolutionary status of V921 Sco from our 
constraints on the disk gas velocity field. The decretion 
disks around post- main-sequence supergiant B[e] stars 
should exhibit a strong ou tflowing velocity component 
(|Lamers fe Pauldrachl I1991D , which is not observed in 
our data. On the other hand, our observation of a Keple- 
rian rotation profile is consistent wit h the expected veloc- 
ity fi eld in a viscous accretion disk (jShakura fe Sunyaevl 
Il973f ) around a pre- main-sequence Herbig B[e] star. 

5.3. On the origin of the Hrj- emission in Herbig Be 

stars 

Arguably, the most important advancement provided 
by our observational data set on V921Sco is that it al- 
lows us to constrain simultaneously the geometry of the 
circumstellar dust disk and the spatial distribution and 
kinematics of the hot gas in the system (see Fig. [7]). Such 
constraints are important, in particular as earlier spa- 
tially unresolved studies have associated the origin of the 
Br7-line emission with fundamentally different physical 
mechanisms, including mass infall (e.g. magnetospheric 
accretion, gaseous inner disks) and mass outflow pro- 
cesses (e.g. stellar winds, X- winds, or disk winds). For 
V921 Sco, we find that the Br7-cmitting gas clearly shows 
a rotation-dominated velocity field and orbits the star 
in the same plane (ip = 145.0 ± 1.3°, i — 54 ± 2°, see 
Tab. 2) as the circumstellar dust disk Up = 145.0 ± 8.4°, 
i = 4 8.8 ± 4°), as i l lustra ted in Fig. [TT] On the other 
hand, iBenistv et al.1 (|2010l) found for the Herbig Be-star 
component in the Z CMa binary system that the By- 
line emission was displaced along the outflow axis, i.e. 
presumably perpendicular to the disk plane. Given that 
the Z CMa observations were obtained during an ex- 
ceptionally strong outbursting phase of the Z CMa sys- 
tem in 2008, they likely do not represent typical condi- 
tions in a Herbig Be star. Another examp le is provided 
by MWC 297, where IWeigelt et ail (|2011[ ) could repro- 
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Fig. 11. — Zoom-in on V921 Sco, covering spatial scales over more than five orders of magnitude. On scales of a few arcminutes (A), 
V921Sco is embedded in an extended ridge-like structure, as shown in this Magellan/IMACS R/[S Il]/H« image (Paper I). On smaller 
scales (B), V921Sco is surrounded by a bipolar nebula with intruiging substructure (C). The axis of the nebula is perpendicular to the 
plane of the disk discovered on scales of a few milliarcsecond within our interferometric image (D). The black ellipse illustrates the location 
of the inner continuum disk truncation radius R ln , as determined with model fitting (Tab. [l] model DISK). Most of the Br7-linc emission 
emerges from inside of the dust sublimation radius, as determined with our kinematical model (E), which is constrained simultaneously by 
the VLTI/AMBER visibilities, VLTI/AMBER+VLT/CRIRES differential phases, and the high-resolution line profile. The image shown 
here is a color-composite of our model channel maps for gas velocities of -35 (blue), (green), and +35 kms" 1 (red). The gas rotation is 
also revealed in a model-independent fashion by the measured photocenter offsets (F). 



duce important features of the measured wavelength- 
differential phases using a magnctosphcric disk wind 
model. These examples, V921 Sco (where Br7 was found 
to trace orbiting gas in the disk plane) , Z CMa (where the 
Br7-emitting gas was observed in a bipolar wind), and 
MWC 297 (where the Br7-emitting might trace a disk 
wind), illustrate that the Br7 emission is probably not as- 
sociated with a single astrophysical process, but can trace 
both mass accretion and outflow processes, as we have 
already concluded from our small VLTI/A MBER Br7- 
survey with spectral resolution R = 1500 (|Kraus et al.l 
l2008bl ). Spatially and spectrally resolved observations 
on a larger object sample will be essential in order to 
explore under which circumstances the Br7-line traces a 
certain mechanism and to determine the relation with 
spectroscopic diagnostics. 

6. CONCLUSIONS 

In this paper, we have investigated the milliarcsecond- 
scale environment around the B[e] star V921 Sco and con- 
strained the spatial distribution and kinematics of ion- 
ized hydrogen gas in the system with a spectral resolu- 
tion up to R = 100,000. We summarize our findings as 
follows: 

• Our H- and A'-band continuum interferometric im- 
ages reveal a spatially extended (Gaussian FWHM 



7.5 mas) disk-like structure, seen under an inter- 
mediate inclination angle of <~ 48.8 ± 4°. Using 
a temperature-gradient model, we determine the 
dust sublimation radius to 1.59±0.25 mas and find 
that the apparent disk size increases with wave- 
length, consistent with an irradiated dust disk. 

• Our VLTI/AMBER (R = 12, 000) and 
VLT/CRIRES (R = 100,000) observations 
spatially and spectrally resolve the Br7-line emis- 
sion from V921 Sco on sub-milliarcsecond scales. 
Using a model-independent photocenter analysis 
technique and our detailed kinematical modeling, 
we find that the line-emitting gas is located in a 
Keplerian-rotating disk, which extends down to 
a few stellar radii. We interpret this finding as 
strong evidence for the pre-main-sequence nature 
of the object, since the decretion disks in post- 
main-sequence B[e]-stars are believed to exhibit a 
significant outflowing velocity component. 

• From our disk kinematical modeling, we derive a 
mass of 5.4 ± 0.4 M Q • (d/1150 pc) for the central 
ob ject. Assuming the distance of 1.15 kpc proposed 
bv lBorges Fernandes et ail (|2007t) , we find that this 
mass is too low to be consistent with the early B- 
type spectral classification, which might indicate 
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that the distance to V921 Sco is probably consid- 
erably larger (~ 2 kpc). Alternatively, the disk 
might be associated with the less massive compo- 
nent in the system. In this dynamical scenario, 
the material from the circumbinary disk is accreted 
onto the intermediate-mass (late B-type) compo- 
nent, preventing significant accretion onto the early 
B-type star. In any case, a careful re-evaluation of 
the spectral classification of the components in this 
enigmatic binary system will be required in future 
studies. 



Our FIRE observations reveal a rich near-infrared 
emission line spectrum of permitted and forbidden 
emission lines (Fe II, [Fe II], C I, He I, O I, N I, 
Mg II, Al II), but show no sign of CO bandhcad 
emiss ion, which s upports a pre-main-sequence na- 
ture (Kraus 2009). 



• From the FIRE spectra, we derive the line flux for 
61 hydrogen recombination lines, including transi- 
tions from the Paschen, Brackett, and Pfund series. 
Utilizing the spatial information provided by our 
VLTI/ AMBER spectro-interferometric data allows 
us to model the line decrements with a hydrogen 
excitation model and to derive electron densities 
for the line-emitting gas. The derived high num- 
ber densities (N e = 2. ..6 x 10 19 m -3 ) suggest that 
the disk around V921 Sco i s particularly massive , 
reinforcing the conclusion of lHenning et al.l (|1998l ) , 
who derived a total gas mass of 40 Mq for the sur- 



rounding millimeter core. 

Our combined spectro-interferometry & line decrement 
modeling approach provides a novel tool to estimate the 
gas density in the innermost AU of protoplanetary disks. 
Applying this tool on a larger sample of T Tauri, Her- 
big Ae/Be, and transitional disks might allow, for in- 
stance, to quantify gas depletion as function of disk evo- 
lution, or to confirm the expected relations between the 
mass accretion rate, dis k surface density, and stellar mass 
(e.g. iCalvet et al.1l200l . 

Furthermore, our study demonstrates the power of 
combining infrared spectro-interferometry with the tech- 
nique of spectro-astrometry. Spectro-interferometry is 
indispensable to characterize the continuum geometry 
(measurement of the dust disk geometry and, for in- 
stance, the detection of companions) and to resolve the 
detailed geometry and kinematics of the circumstellar 
gas. On the other hand, spectro-astrometry provides a 
ressourcc-cfficient and straightforward method to mea- 
sure first-order kinematical information with very high 
spectral dispersion in the spatially unresolved regime, 
providing highly complementary constraints for the in- 
terpretation of spectro-interferometric signatures. 
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